Weak Interaction Rate Coulomb Corrections in Big Bang Nucleosynthesis 
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We have applied a fully relativistic Coulomb wave correction to the weak reactions in the full 
Kawano/Wagoner Big Bang Nucleosynthesis (BBN) code. We have also added the zero temperature 
radiative correction. We find that using this higher accuracy Coulomb correction results in good 
agreement with previous work, giving only a modest ~ 0.04% increase in helium mass fraction over 
correction prescriptions applied previously in BBN calculations. We have calculated the effect of 
these corrections on other light element abundance yields in BBN and we have studied these yields 
as functions of electron neutrino lepton number. This has allowed insights into the role of the 
Coulomb correction in the setting of the neutron-to-proton ratio during the BBN epoch. We find 
that the lepton capture processes' contributions to this ratio are only second order in the Coulomb 
correction. 

PACS numbers: 14.60.Pq; 14.60.St; 26.35.+C; 95.30.-k 
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I. INTRODUCTION 

The study of Big Bang Nucleosynthesis (BBN) has 
been and is a powerful tool for testing cosmological mod- 
els and constraining the fundamental parameters of the 
universe. Since primordial nucleosynthesis occurs rela- 
tively soon after the Big Bang (~ Is), BBN provides one 
of the best windows into the physics of the early universe. 

Before the launch of the Wilkinson Microwave 
Anisotropy Probe (WMAP), BBN predictions along with 
direct observation of the primordial element abundances 
were used to constrain the baryon-to-photon ratio, 77. 
The independent high precision determination of 77 from 
the ratio of the acoustic peak amplitudes in the cosmic 
microwave background (CMB) from WMAP [H, i, Q al- 
lows us now to use BBN to constrain other unknowns 
in the early universe and physics beyond the Standard 
Model. 

The CMB measurement of 77 increases in precision with 
accumulating WMAP data and the future Planck mission 
promises even higher precision, with a projected ~ 1% 
accuracy in ry [J 0] • The measurement of primordial 
deuterium also shows promise for higher accuracy deter- 
mination as more QSO lines of sight become available 
[^, Q . This fuels the motivation to further refine the 
calculation of predicted primordial element abundances. 
For this reason we have analyzed the effect of adding the 
full relativistic Coulomb wave correction factor (relativis- 
tic Coulomb barrier factor) to the weak reaction rates in 
the BBN calculation. 

Among many issues, a key piece of physics that sets 
the stage for primordial element nucleosynthesis is the 
evolution of the neutron-to-proton ratio, n/p. The n/p 
ratio is critical in determining the synthesis of the pri- 
mordial elements because it sets the number of neutrons 
available to build nuclei. 

The neutron-to-proton ratio is effectively determined 
by the competition between the charge-changing weak 
interaction rates and the expansion rate of the universe. 
Listed below are the weak reactions which interconvert 



neutrons and protons: 
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The corresponding rates for these weak reactions are 
denoted by A^^„ and Ae-p, Ap^p and Ae+„, An^^^^^ and 
Ape-pe for the forward and reverse reactions in Eq. ([1]), 
Eq. ([2]), and Eq. ([3]), respectively. Defining A„ = Ap^„ -I- 
Ae+n + '^iidccay ^-iid Atot = A„ -f Ap^ -I- Ae-p + Apg-p^, and 
defining the neutron-to- proton ratio to be n/p, we can 
show that in the early universe 
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where t is the Freidmann-Lemaitre-Robertson- Walker 
timelike coordinate 0]. Note that the net number of 
electrons minus positrons per baryon is Yg = (n-e- — 
rie+)/nf, = (1 -I- n/p)~^. At high temperatures, T » 1 
MeV, the weak reaction rates are fast compared to the 
expansion rate of the universe, steady state equilibrium 
{-^{n/p) = 0) is a good approximation, and the neutron- 
to-proton ratio is given by 0, HI 



n ^l^eP Ag-p + Apg-^^ 

P Ap^n Ag+^^ -\- An decay 



(5) 



The relatively slow rates at high temperatures for both 
free neutron decay and the corresponding reverse three 
body reaction allow the steady state equilibrium neutron- 
to-proton ratio to be approximated as 
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If the neutrinos have thermal, Fermi-Dirac energy distri- 
bution functions, this can be written by 



n _ (Ae-p/Ae+„) +6-'^-.+"-"^ 

P ~ (Ag-p/Ag+nje^-e-^'+f + l' 
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FIG. 1: Neutron-to-proton ratio as a function of temperature. 
The full standard BBN zero lepton number case is given by 
the solid line. The dashed line is the neutron-to-proton ratio 
as calculated with an enforced assumption of steady state 
equilibrium, i.e., Eq. (|8|l. 



where r],y^ = Hu^ /T is the electron neutrino degener- 
acy parameter, rj^ — fi^/T is the electron degeneracy 
parameter, with and /ie the electron neutrino and 
electron chemical potentials, respectively, and ^ is the 
neutron-proton mass difference divided by temperature, 
^ = Sninp/T with Srrinp = to„c^ — nipC^ w 1.293 MeV, 
where we take the Boltzmann constant to be /cb = 1 01 ■ 
At high temperatures the weak reactions can be fast 
enough to maintain chemical equilibrium. In chemi- 
cal equilibrium, the chemical potentials satisfy the Saha 
equation, /i^^ -\- Hn = He- + Mp- When chemical equilib- 
rium is maintained the neutron-to-proton ratio will be 



P ^ 



(8) 



This result can be obtained directly from the ratio of the 
appropriate Fermi-Dirac distribution functions or, alter- 
natively and equivalently so long as all reactants have a 
Fermi-Dirac form for their energy spectra, from evalua- 
tion of the rates in Eq. Q 0]. 

As the universe cools and expands, the rates of the 
weak reactions become slow compared to the expansion 
rate of the universe. At this point chemical equilibrium 
can no longer be maintained and a period known as the 



"weak freeze out" occurs as the neutron-to-proton ratio 
pulls away from its equilibrium value. For a faster expan- 
sion rate of the universe the weak reaction rates become 
comparatively slow earlier and the neutron-to-proton ra- 
tio falls out of equilibrium at higher temperatures, yield- 
ing a relatively larger n/p value. Likewise, for a slower 
expansion rate, the weak reactions can maintain chem- 
ical equilibrium longer and the neutron-to-proton ratio 
consequently would be lower. The expansion rate of the 
universe is set by the local total energy density through 
the Freidman equation. 

Figure[T] shows the neutron-to-proton ratio as a func- 
tion of temperature for a standard Big Bang scenario 
with zero lepton numbers, i.e., /i^^ — /ip^ — fie ^ 0- This 
figure shows both the actual n/p ratio and the approx- 
imation to this with an enforced chemical equilibrium 
condition. Obviously, these agree for high temperature 
but diverge once the weak reaction rates become slow 
compared to the expansion rate of the universe. Note 
that the actual n/p ratio becomes constant once nearly 
all free neutrons are incorporated into alpha particles at 
T < 100 keV. This figure shows the n/p ratio only for 
free neutrons and neutrons bound within alpha particles, 
and neglects the neutrons bound in ^H, ''H, •^He, and nu- 
clei heavier than ^He, accounting for the small divit in 
the upper curve near T =80 keV. 

Primordial element abundance yields are calculated by 
a BBN code that time evolves the temperature and ex- 
pansion rate of the universe along with the nuclear and 
weak reactions rates. We have used a modified version 
of the Kawano/ Wagoner BBN code @, El [IH, [ll in or- 
der investigate the effect of integrating the relativistic 
Coulomb barrier factor in the appropriate weak reaction 
rates. In section II we discuss the calculation of the weak 
reaction rates. In section III we discuss the relativistic 
Coulomb correction employed here and Coulomb correc- 
tion prescriptions studied previously. In section IV we 
present results and give a discussion and in section V we 
give conclusions. 



II. THE WEAK REACTION RATES 

We calculate the individual weak interaction rates with 
the following phase space factor forms and with a com- 
mon matrix element which is proportional to the inverse 



of an effective ft-value (ft) 
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In 2 



/ {Ee + 5mnpfEe(El - mec^f'^ [S^+] [1 - 5pJ dEe, 



(11) 



In 2 



/ F [Z, Ee] {Ee - 6mnpfEe{E^, - mecY^^ [S.J [1 - S,-]dEe, 



(12) 



In 2 



F [Z, Ee] {Smnp - EefEe {E^ - TTleC^) [1 - 5pJ [1 - Se-]dEe, (13) 



^ _ In 2 r^"""^ 



F [Z, Ee] (Smnp - EefEe {E^ - UleC^Y'^ [5pJ [S^-] dEe 



(14) 



where Ee is the total electron or positron energy as ap- 
propriate, TrieC^ is the electron rest mass, and F [Z, Ee] 
is the Coulomb correction Fermi factor which will be dis- 
cussed in detail below. Note that the nuclear charge rele- 
vant here is Z = 1. Se-/+ and S'l/^/p^ are the phase space 
occupation probabilities for electrons/positrons and neu- 
trinos/antineutrinos, respectively. For neutrinos and 
electrons with energy distributions with the expected 
thermal form, the occupation probabilities are 
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where Ti, is the neutrino temperature parameter, r/i, is 
the neutrino degeneracy parameter (the ratio of chemi- 
cal potential to temperature), and E^ is the appropriate 
neutrino or antineutrino energy. 
We take 



In 2 _ c(mec^)' p (l + 3|CA/Cy| 

(ft) ~ he ' ' 27r3 



(18) 



where Gp ~ 1.166 x 10~^^MeV~^ is the Fermi con- 
stant, Cv and Ca are the vector and axial vector cou- 
pling constants, respectively, and we have taken the ab- 
solute squares of the Fermi and Gamow- Teller matrix 
elements for the free nucleons to be |Afi?p = 1 and 
\Mgt\'^ = 3, respectively. Here 5 is a factor which in- 
cludes both Coulomb and other ( "radiative correction" ) 
effects which amount to a few percent change in the ef- 
fective /t- value, (ft). 

Of course, Cy and Ca are coupling constants that are 
renormalized by the particular strong interaction envi- 
ronment characterizing free neutrons and protons. (Ab- 
sent strong interactions Cy = Ca = 1-) Given that these 



are a priori unknowns, as is (5, we follow the standard pro- 
cedure 12]: we take the free neutron decay rate as the 
product of Eq. pB]) and the phase space factor in Eq. 
(with Si}^ — Se- — 0) and we then set this equal to the 
inverse of the laboratory-measured free neutron lifetime, 
T„. The world- average of the laboratory measurements 
is r„ = 887.7 seconds [Tsj . 

Note that changing the prescription for the Coulomb 
correction factor F[Z, Eg] in Eq. (|13p will have the ef- 
fect of renormalizing the effective free nucleon weak in- 
teraction matrix elements {i.e., renormalizing (ft)) for a 
given r„. As we will see below, this renormalization will 
be the dominant component of the Coulomb correction 
alteration in, e.g., the ^Hc BBN yield. 

The rates for all the individual weak reactions are 
shown as functions of temperature in Fig. O At high 
temperatures the forward and reverse rates of the lep- 
ton capture reactions in Eq. ([T]) and Eq. ^ dominate 
the neutron- proton inter-conversion process. Note that 
the rates for the forward process in Eq. ^ and the re- 
verse process in Eq. ([TJ are affected by the threshold, 
Smnp+meC^ . At lower temperatures this threshold makes 
these rates relatively slower than the rates for the lepton 
capture channels without this threshold, i.e., the forward 
process in Eq. ([T]) and the reverse process in Eq. ([2]). 

This figure shows that at a lower temperature (T ^ 
SrUnp) the electron capture rate Ag-p and the three-body 
rate \pe-u^ track each other closely, differing by a factor 
of order unity. This is readily explained as follows. First 
note that the integrands in the phase space factors in 
Eq. ^ and Eq. (|14p are identical. Although the former 
phase space factor is proportional to 1 — S,y^ and the latter 
to Sjj^ , when the Ve and De energies in these distributions 
are expressed in terms of electron energy, E^^ — Ee — 
Sninp and Ep^ — Srrinp — Ee, respectively, we see that 
1 — iSiy^ — Sjj^ . Second, though the limits of integration 
for these phase space factors differ, we note that when 
T <C SrUnp the upper limits are effectively the same. Only 
the lower limit is different in the two cases, « 1.3 MeV in 
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FIG. 2: All six weak reaction rates as a function of tempera- 
ture. The solid (red) line is for \v^n, the dashed (green) line 
is for A(.+„, the dotted (blue) line is for An^^^^y, the small- 
dashed(pink) line is for \pv^, the dotted-dashed (cyan) line is 
for Ag-p, and the black dotted-spaced line is for ^^p^-n^- All 
lepton chemical potentials are set to zero here. 



the former and w 0.511 MeV in the latter. 

As the temperature decreases, free neutron decay be- 
comes the dominant weak reaction. This remains the 
case through the epoch when strong and electromagnetic 
nuclear reactions freeze out of equilibrium ("nucleosyn- 
thesis"), when nearly all free neutrons are incorporated 
into alpha particles. 



III. COULOMB CORRECTION TO THE WEAK 
INTERACTION RATES 



For weak interaction processes that have an electron 
and a proton in either the initial or final state, the 
Coulomb interaction must be taken into account. In fact, 
the phase space factors presented above are derived by us- 
ing plane wave functions for the entrance and exit chan- 
nel leptons, but then "correcting" where Coulomb waves 
should be used instead by multiplying the appropriate 
phase space integrals by the Fermi factor F{Z, Eg). 

The Coulomb potential is attractive in the e^/p chan- 
nel. This has the effect of increasing the electron prob- 
ability amplitude at the nucleus (proton) and, in turn, 
this will always increase the affected phase space factors. 
In other words, the phase space factors in the expressions 
for the rates for both the forward and reverse processes 
in Eq. ([1]) and Eq. ^ will be increased over a case where 
only plane waves are used. 



A. Previous Corrections to the BBN Weak 
Reaction Rates 

The Coulomb correction to the weak rates in BBN was 
first employed by Wagoner in an early version of 
the BBN code. Wagoner took a representative value of 
the correction from around the time of weak freeze out, 
when the weak rates have the largest effect on the n/p 
ratio, and used this to scale the neutron lifetime, t„. 
This had the effect of increasing the effective neutron 
lifetime over what had been chosen for t„ by about 2%. 
This effectively slowed down all weak interaction rates 
by 2% because the weak rates are all normalized by the 
neutron lifetime. Wagoner's correction over-estimated 
the Fermi factor, giving an erroneous increase in helium 
mass-fraction, Yp, of about 0.5%. This over production 
was largely a result of "correcting" weak reactions that 
should not have been corrected, e.g.^ n + ^ p + P,,- 
Wagoner did not include any radiative corrections. 

Dicus et. al [l^ were the first to add an energy 
dependent Coulomb correction to the BBN calcula- 
tion, along with the zero-temperature radiative correc- 
tions and finite-temperature radiative corrections. They 
approximated the Coulomb correction using the non- 
relativistic form for the Fermi factor 
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(19) 



where (3 — w/c is the electron velocity. It was pointed 
out in Ref. ^] that Ref. p"], like Wagoner, Coulomb- 
corrected the rates that should not have had a Fermi 
factor in their phase space integrands. 

Ref. [isj derived the zero-temperature radiative correc- 
tions for a point nucleon, finding that all the weak rate 
integrands should be multiplied by 
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where (3 is again the electron velocity and y and e are 
the neutrino energy and electron energy divided by the 
electron mass, respectively, and 

C{(3,v) = 40-)-4(C/-l)(y/3e-3/2-|-ln2y) (21) 
+U {2{l + (3^) + / Qt^ -'i(3U\ 
-4[2 -I- 11/3 -I- 25/3^ + 25/3^ + 30/3'' 
+20/35 + 8/36)/(H-/3)^ 



Here U is defined to be 



U EE /3"4anh~^/3. 



(22) 



The corrections in Ref. jl^] resulted in a ~ 0.4% reduc- 
tion in Yp from a calculation with the Coulomb effect 
alone, plus a ~ 0.2% increase stemming from the zero- 
temperature radiative corrections. 

Ref. [Ill and Ref. [1^ appropriately applied the 
Coulomb correction to only those rates which require 
one. Ref. [2l| used the Fermi factor approximated at 



5 



order a, while Ref. used a non-relativistic version of 
the Fermi factor in Eq. p^ . They also applied the zero- 
temperature radiative corrections defined above as well 
as several other corrections. 



B. New Coulomb Correction and Modifications to 
the BBN code 



Additionally, this is the first time this version of the 
correction has been applied to the full reaction network 
in the Kawano/ Wagoner code. This allowed us to study 
the effect of the correction on all of the light element 
abundances. We were also able to study the effect of the 
correction on nucleosynthesis in the presence of neutrino 
degeneracy (a lepton number). 



In this work we employ a version of the Coulomb cor- 
rection which can better take into account the poten- 
tially relativistic kinematics of initial or final state elec- 
trons. We use the Coulomb correction that is discussed 
in Ref.s fTllllIll: 



G{±Z,Ee)=xF{±Z,Ee), 



(23) 



and we define x = [E^ — {meC^YY^'^ / Eg, the ratio of 
charged lepton (electron here) momentum to energy. In 
Eq. P5|) . F{±Z,Ee) is the Fermi factor (or relativistic 
Coulomb barrier factor) approximated here as 



F{±Z,E,) ^ 2(1 + s)i2pR) 



2(s-l)g7n^ 



r(s -I- iuj) 



Ti2s + 1) 



(24) 



In this expression, the upper signs are for electron emis- 
sion and capture, the lower signs are for positron emission 
and capture in the general case for a nucleus of electric 
charge Z (in our case Z = 1), s = [1 — {aZ)^]^^^, a is the 
fine structure constant, uj = ±Z/x ("-f" for the in our 
cases), and R is the nuclear radius in electron Compton 
wavelengths, R = 2.908 x lO'^A^/^ - 2.437^-^/^ where 
A is the nuclear mass number and A=l in our case. This 
is the most accurate Coulomb correction that has been 
employed in a BBN calculation. 

In order to properly apply these features of the cor- 
rection, we used a version of the Kawano/ Wagoner code 
where the weak reaction rates have been entirely rewrit- 
ten. We will only briefly describe this code here. A 
detailed description can be found in Ref. [sj. 

The original Kawano/ Wagoner code calculates the 
weak rates with a total lumped sum of the n — > p and 
p — > n rates: 

An = ^}y^+n^p+e- + Ki+e+^p+P^ + ^n^p+e-+P^ (25) 



Xp Ap_(_e- 



A 



(26) 



In our version, we have separated these summed rates to 
calculate all 6 weak reaction rates individually. Another 
key feature of the code is that each rate is modularized, so 
that any neutrino and antineutrino distribution function 
and time dependence thereof can be applied. 

Because of this modularization we were able to apply 
a Coulomb correction to only those rates that require 
one. In other words, we were also able to include an 
appropriate relativistic Fermi factor in the integrand of 
those weak rates. 



IV. RESULTS AND DISCUSSION 

We have applied the Coulomb correction described 
above along with zero temperature radiative corrections 
in the full Kawano/ Wagoner BBN code. The integrated 
effect of these corrections can be seen by the changes in 
the light element abundances. 

The Coulomb correction described above affects BBN 
abundance yields in a subtle, but interesting way which 
gives insight into the weak interaction's role in setting 
the neutron abundance in the early universe. First, as 
outlined in the last sections, the key effect of calculat- 
ing the weak rates with Coulomb waves instead of plane 
waves is to increase the electron's probability amplitude 
at the proton. This means that the rates corresponding 
to e~ -f p — > n -I- j/g, Ve + n ^ p + e~ , n ^ p + e~ + v^, 
and p + + De ^ n will all increase over plane-wave 
calculated rates. 

This is true, but in the BBN calculation the net ef- 
fect of adding a Coulomb or radiative correction which 
increases the phases space factor for free, vacuum neu- 
tron decay is to reduce the weak matrix element [increase 
{ft)) common to all of the rates of the forward and re- 
verse processes in Eq. ([T]), Eq. and Eq. This is 
because for a given vacuum {Se = S^^ = 0) neutron life- 
time, r„, we set An^^^^^ | vacuum = T~\ and an increased 
phase space factor then implies an increased value for 
(/<)• 

Therefore, the chief effect of a Coulomb correction- 
mediated increase in phase space factors is a decrease in 
the overall strength of the weak interaction. In turn, a 
weaker weak interaction would cause a higher tempera- 
ture for freeze-out from chemical equilibrium and a con- 
comitant increase in the neutron-to-proton ratio emerg- 
ing from the weak freeze out process. Since, to first ap- 
proximation, all neutrons will eventually be incorporated 
into alpha particles, the phase space factor-enhancing 
Coulomb correction should give rise to an increase in the 
primordial mass fraction Yp. 

In broad brush this is indeed what the BBN calcula- 
tions show. Our modularized code, which allows us to 
track the individual weak rates, affords us a deeper in- 
sight into what is happening. Though all rates are renor- 
malized downward by the Coulomb correction, the rates 
for the particular processes with a Fermi factor in their 
phase space integrals (Eq.s (H)), HI]), (US]), and HI])) are 
decreased less. In other words, they are increased rela- 
tive to the rates in Eq. PH]) and Eq. pT]) . Ap^p and Ae+„, 
respectively. 
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TABLE I: ''He mass fraction Yp and deuterium abundance 
D/H as calculated with our code for various implementations 
of Coulomb and radiative corrections as indicated. The Base- 
hne table entries are the uncorrected values; table entries des- 
ignated by Wagoner, Esposito, and Lopez and Turner were 
com put ed using the correction prescriptions in Ref. [2H . 
and [23I , respectively, but with our code and with the current 
world-average neutron lifetime. 





Yp 


D/H 




Baseline 


0.239 


2.522 X 10"^ 




Wagoner 


0.2427 


2.543 X 10"^ 


1.16% 


Esposito 


0.2416 


2.537 X 10-^ 


1.09% 


Lopez and Turner 


0.2416 


2.537 X 10"^ 


1.09% 


New Correction 


0.2417 


2.543 X 10~^ 


1.13% 


New Correction 


0.2422 


2.542 X 10"^ 


1.34% 


with Zero- Temperature 








Radiative Corrections 









Nevertheless, neutron decay has more leverage over the 
eventual n/p ratio than do the lepton capture processes. 
At higher temperatures where the n/p ratio is well ap- 
proximated as the ratio of the sum of the neutron pro- 
duction rates to the sum of the neutron destruction rates, 
Eq. ([HI), note that the small and comparable fractional 
relative increases in Xe-p and Xv^n tend to compensate 
each other to first order. This is because Xf,-p is in the 
numerator and X^^n is in the denominator in Eq. 
As a consequence, changes in n/p stemming from lep- 
ton capture processes are second order in the Coulomb 
corrections. 

Table I presents a comparison of BBN calculations 
of the *He mass fraction, Yp, and the deuterium abun- 
dance relative to hydrogen, D/H, all performed with our 
code. Shown in this table are a "Baseline," standard 
(/iy^ = /ie = 0) BBN case with no Coulomb corrections 
(i.e., F{Z,Ee) = 1), and cases where the same calcula- 
tions were done but where the Coulomb corrections of 
Wagoner Ref. fiol. E sposito et al. Ref. [2l|, and Lopez 
and Turner Ref. [2^, were used. This table also shows 
results from the same BBN calculation but using our 
new relativistic Coulomb correction for cases with and 
without a radiative correction. Consistent with the ar- 
guments given above, we see that the Coulomb correction 
prescriptions in all of these cases give a ~ 1% increase 
in Yp. All of the various versions of the Coulomb cor- 
rection are consistent with each other. The relativistic 
Fermi factor used in this work results in a modest 0.04% 
increase in Yp over the result in Ref. [2l| . 



The modular nature of the weak rates in our BBN code 
allows us to examine the effects of the Coulomb correc- 
tion for scenarios in which the lepton numbers residing 
in the v^. and seas are nonzero. This is the first such 
study of the Coulomb and radiative correction effects in 
a case with nonzero values for ^.^^ and /.tp^ . (Here we 
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FIG. 3: Percent change in nuclear abundances from the Base- 
line values as a function of electron lepton number from when 
the Coulomb and zero-temperature radiative corrections are 
included. 



take + /Xp^ = 0, reflecting assumed neutrino chemical 
equilibrium at high temperatures.) We define the lepton 
number residing in these neutrino seas to be 



(27) 



and 



Ve proper num- 



IS 



where n^^ and rip^ are the v, 
ber densities, respectively, and — (2^(3)/7r 
the corresponding photon number density with ^(3) k, 
1.20206 the Riemann-Zeta function of argument 3. 
The primordial helium abundance plus observationally- 
and experimentally-determined neutrino flavor oscilla- 
tion data restrict |L^J < 0.1 i26i, |23, IH . This up- 
per limit conceivably could rise to ~ 0.2 if allowance is 
made for additional contributions to the energy density 
in the early universe [1^, [s^]- Models which attempt 
to reconcile light-mass sterile neutrinos with BBN and 
large scale structure plus cosmic microwave background- 
derived overall neutrino mass closure constraints usually 
invoke lepton numbers. But the lepton numbers invoked 
in these schemes can increase the *He-based upper limit 
on \ Lu, I 3l [ . 33. Sterile neutrino dark matter scenarios 



[il HIMliM il, S ii, H Si E ^ 



|3a, iSg, IM |4aM, IM lia I46|, |47 

also can invoke significant lepton numbers. We therefore 
consider a range < L^^ < 0.3 as an interesting example. 

Unlike the neutrino degeneracy parameter, Li,^ is not 
a co-moving invariant. The ratio of the neutrino temper- 
ature parameter, T,y, in and S'p^ to the temperature 
of the plasma, T^, evolves in time. This ratio starts out 
as unity for T 1 MeV and, as electrons and positrons 
annihilate and transfer their entropy preferentially to the 
photons and plasma, eventually asymptotes to (4/11)^^'^ 
at low temperatures. The lepton numbers and the neu- 
trino degeneracy parameters are related by 



(28) 
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and at small lepton number rj^^ « lA&L^^. Here we 
consider only positive values of L^^, i.e., cases with a 
preponderance of UeS over v^s, as these are the most 
interesting with respect to '*He. 

An effect of nonzero /i^^ and /ip^ will be to change the 
energy weighting in the phase space integrands in Eq. ([9]), 
Eq. inD, Eq. ((T3]), and Eq. We might then expect 

a concomitant alteration in the effect of the Coulomb 
correction over the zero lepton number cases for Ag-p, 
Ai/^n, An^^^^y , and Apg-p^, respectively. In Fig. [3] we show 
the relative change in BBN abundance yields for ^He, 
^Li, ^H, ^H, and "^He, over the Baseline (no Coulomb 
or radiative corrections) case as a function of L^^. In 
the figure the percent change in the '*He mass fraction 
is designated by 5*IIe, while the percent change in the 
abundances for deuterium, tritium, and ^He are given by 
(5D/H, (5T/H, and (5^He/H, respectively. First, we note 
that the overall sense of the Coulomb correction is to 
increase Yp and all of the nuclear abundances yields over 
the Baseline case for the lepton number range examined. 

A higher and the accompanying higher number 
density Uij^ and lower number density rtp^ , will have the 
effect of increasing Xj^^n and A„j^^^y and decreasing Ae-p 
and Ape-pj over the zero lepton number case. This is sim- 
ply a result of an enhancement or reduction in final state 
blocking or entrance channel lepton number density as 
appropriate. Figure 3 shows that while the '^He and '^He 
Coulomb correction abundance yield enhancements are 
essentially flat with increasing L^^ , the abundance yield 
enhancement for ^Li increases while the corresponding 
enhancements for '^H and decrease with increasing 
L^^. The trends with L^^ of ^Li/H, ^H/H, and ^H/H, 
versus those for the ^■"'He yields reflect the different times 
at which these species are formed and the sensitivity of 
the relevant reaction production mechanisms to the local 
neutron abundance and temperature. 



work of the Kawano/ Wagoner BBN code. We have used 
this code to study BBN abundance yields for a range of 
neutrino chemical potentials. We find that the fully rel- 
ativistic Coulomb correction essentially agrees with pre- 
vious non-relativistic prescriptions, giving only a 0.04% 
increase in the ''He yield over that in Esposito et al. [2l| . 
Our calculations show interesting trends in the light el- 
ement abundance yields with increasing electron lepton 
number. The modularization of the individual weak in- 
teraction processes in our code has allowed us to gain 
insights into how the rates for these processes are al- 
tered by the Coulomb and radiative corrections and how 
these alterations affect the neutron-to-proton n/p ratio in 
the early universe during the BBN epoch. In particular, 
we point out that the lepton capture processes produce 
changes in the n/p ratio which are only second order in 
the small Coulomb corrections. 

As the accuracy of measurements of the CMB and 
the primordial abundances of the light elements increase, 
BBN will give even better constraints on physics in the 
early universe. Currently, the increase in precision gained 
from including this rclativistic version of the correction 
is probably unnecessary. However, in the future as the 
measurements for the main parameters affecting BBN, 
such as the baryon-to-photon ratio and the neutron life- 
time, increase in precision, it may be beneficial to include 
this version of the correction. 
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